An orbital-selective spin liquid in a frustrated heavy fermion spinel LiV204 
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The pronounced enhancement of the effective mass is the primary phenomenon associated with 
strongly correlated electrons. In the presence of local moments, the large effective mass is thought to arise 
from Kondo coupling, the interaction between itinerant and localised electrons. However, in d electron 
systems, the origin is not clear because of the competing Hund's rule coupling. Here we experimentally 
address the microscopic origin for the heaviest d fermion in a vanadium spinel UV2O4 having geomet- 
rical frustration. Utilising orbital-selective 51 V NMR, we elucidate the orbital-dependent local moment 
that exhibits no long-range magnetic order despite persistent antiferromagnetic correlations. A frus- 
trated spin liquid, Hund-coupled to itinerant electrons, has a crucial role in forming heavy fermions with 
large residual entropy. Our method is important for the microscopic observation of the orbital-selective 
localisation in a wide range of materials including iron pnictides, cobaltates, manganites, and ruthnates. 

PACS numbers: 



Introduction 

Electrons in metals behave as quasiparticle dressing inter- 
actions. The mass of quasiparticles often becomes extremely 
heavy when metallic phases are close to the quantum critical 
boundary for the insulating or magnetic phase 1 . The micro- 
scopic understanding of heavy quasiparticles (HQs) has been 
a goal of modern many-body quantum statistics. An estab- 
lished route for HQs is antiferromagnetic Kondo coupling be- 
tween localised / spins and itinerant electrons in rare-earth 
metals. In contrast to the /-electron case, the presence of lo- 
calised spins is not apparent for d-electron systems. Alter- 
native routes driving d HQ formation have been challenging 
issues in strongly correlated electron physics. 

A representative d HQ material is the vanadium spinel 
L1V2O4 (Ref. 2. 3), which has a highly-frustrated pyrochlore 
lattice for the B site V 3 5+ (3d 15 ) ions (Fig. la). Anisotropic 
orbital-dependent intersite interactions give an itinerant e' g or- 
bital and a more localised a\ g orbital through a small trigo- 
nal distortion of the VC>6 octahedron (Fig. lb) 4,5 . The HQ 
was initially explained by off-site Kondo exchange interac- 
tions, /k, between localised a\ g moments and itinerant e' g 
electrons (Fig. lb)i. In the tig manifold, however, the strong 
on-site ferromagnetic Hund's exchange interaction, Jn, can 
overcome J^. Many alternative scenarios such as geometrical 
frustration via antiferromagnetic interactions Jaf, 6-8 , electron 
correlations^*^, and spin-orbital fluctuation s 11 ! 12 have been 
proposed. 

Experimentally, the interpretations have been unclear 
for HQs in L1V2O4. Charge- sensitive probes such as 
resistivity 2 , photoemission 1 — , and optical^ measurements 
showed crossover from a high- temperature incoherent metal 
to a low-temperature Fermi liquid across the characteristic 
temperature r* ~ 20-30 K. In contrast, spin- sensitive probes 
including static spin susceptibility 2,15 and inelastic neutron 
scattering 16 measurements imply local moments with anti- 
ferromagnetic correlations at low temperatures. Furthermore 
anomalous temperature T dependences on the specific heat 
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FIG. 1 : V pyrochlore lattice and exchange interactions between 
3d orbitals in LiV204. (a) Corner- sharing V-tetrahedra offer ge- 
ometric frustration against magnetic ordering. In the cubic lattice 
with local trigonal distortion, three vanadium sites with aig-like or- 
bitals, numbered 1, 2, and 3 are not equivalent to the magnetic field 
Hq rotated from the [001] to [110] axis, (b) A trigonally distorted 
VOg ligand field splits a five-fold 3d level into e g and t2 g having 
lower a\ g and upper e' g with 1.5 electrons. For the 2 and 3 sites, the 
trigonal axis is located at 54° measured from the crystal [001] axis, 
(c) In a Kondo lattice having localised a\ g and itinerant e' g electrons 
with transfers t, exchange interactions between the occupied orbitals 
are composed of the on-site ferromagnetic Hund's coupling Jn, the 
off-diagonal Kondo coupling 7k > and the off- site antiferromagnetic 
kinetic exchange 7af- 



C/T and the Hall coefficient conflict with those expected in 
a conventional Fermi liquid 2 . Despite the theoretical view of 
orbital- selective interactions, no experimental effort has been 
made to detect the orbital degrees of freedom. 

Here, we address the first experimental approach for micro- 
scopic observations of the d HQ via orbital-resolved nuclear 
magnetic resonance (ORNMR) measurements in L1V2O4. 
The previous NMR experiments using closed-shell Li sites^A^ 
only measured the net spin susceptibility proportional to the 
bulk value because the hyperfine interactions at the Li site 
surrounded by 12 vanadium sites average out the anisotropy. 
Our ORNMR spectroscopic approach using on-site 51 V spins 
on a high-quality single crystal is sensitive to the orbital- 
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dependent local spin susceptibility, which is beneficial for 
probing strongly correlated electrons with the orbital degrees 
of freedom. 



Results 

ORNMR Knight shift. In LiV 2 4 the spin susceptibil- 
ity X s consists of the a\ g and e' g components: X s = % a + % e 
(hereafter the superscripts a and e denote a\ g and e g , re- 
spectively). The NMR frequency shift called the Knight 
shift, K s = (K*,Ky,K*) (x, y, and z are the principal axes), 
measures the spin susceptibility via the hyperfine interaction 
Jt^ = Li I " A; • s ; with the hyperfine coupling tensor A; be- 
tween the nuclear spin I and the paramagnetic spin polarisa- 
tion Si for the z-th electron under an external magnetic field. 
Whereas the isotropic shift K? so = (K* +2K%)/3 due to the 
core polarization is proportional to # s , the anisotropic part 
K^ x = 2(K* —K%)/3 due to the orbital- specific dipole hyper- 
fine interaction 17 is expressed by using the principal z compo- 
nent of the coupling constants, A a z and Af , 

Kl = -±-{A a zX a +KX% (1) 

where N is Avogadro's number and /Ib is the Bohr magneton. 
In contrast to K? so , K^ x measures the hyperfine- weighted av- 
erage of the orbital-dependent spin susceptibility. In the ionic 
limit, A z and A e z are given by a quadratic combination of the 
angular momentum (see Methods) with the reversed sign and 
the same amplitude, A a z = —A\ > (Ref. |Tt|). Hence, we 
can distinguish which orbital dominates the spin susceptibil- 
ity from the sign of K^ x and obtain the orbital occupation from 
the amplitude. Namely, K^ x should be positive (negative) for 

x a >x e (x a <x e y 

The 51 V Knight shift tensors of L1V2O4 are determined 
from the angle dependence of K for 2-300 K, as shown in 
Fig. 2. The 51 V NMR spectra were detectable only at specific 
angles, where the nuclear quadrupole interaction almost van- 
ishes, because the nuclear spin- spin relaxation times at other 
angles are too fast. The obtained K traces three cosine curves 
of Eq. (2) in the Methods, which satisfies the cubic Fd3m 
lattice. The principal z axis of K at ±54° for the two V sites 
indicates the 3d orbital symmetry governed by the trigonal 
V0 6 crystal field (Fig. lb). At 300 and 2 K the relationship 
K z > K* (K^ x > 0) shows the ai g -dominant spin susceptibility 
(X a > X e )- The result is consistent with the localised nature of 
the a\ g orbital, as theoretically suggested^^i^. 

Temperature dependence of orbital occupations. To ad- 
dress HQ formation, we measured the thermal variations of 
the 51 V Knight shifts K? so and —K^ x in comparison with the 
7 Li Knight shift 7 K S and the bulk spin susceptibility % (Fig. 
3a). Good linearity was observed between these Knight shifts 
and x (Fig- SI). All of data show a Curie- Weiss-like increase 
at high temperatures, followed by a broad maximum at ap- 
proximately 20 K. The results agree with the spin suscepti- 
bility for high-quality crystals, free from a Curie-tail increase 
at low temperatures 2,3 . The on-site 51 V Knight shift probes 
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FIG. 2: Angle dependence of the 51 V Knight shift K of LiV 2 4 at 
(a) 300 K and (b) 2 K. K is plotted against the magnetic field direc- 
tion measured from the [001] axis to the orthogonal [110] axis. The 
K data with different colors and symbols arise from three vanadium 
sites 1, 2, and 3 due to different orbital configurations against the ex- 
ternal magnetic field and trace three curves of the general formula, 
as shown in Eq. (2). Error bar is defined by a standard deviation of 
the Lorentzian fit of the NMR spectra. The principal Z and X axes of 
K for 2 and 3 are respectively located at =f54° and ±36° measured 
from [001]. 



the spin susceptibility with greater sensitivity than the off- site 
7 Li one and shows a smooth decrease below T*. Below 5 K 
K s becomes nearly T -independent, as observed in the Fermi 
liquid. 

When a\ g local moments becomes Fermi liquid with HQ 
via a\ g -e' g hybridisation or intersite Kondo coupling below T*, 
X a is expect to decrease significantly, whereas x a is l ess sen_ 
sitive. It could lead a decrease of K^ x from Eq. (1). To inspect 
this property, we plot Kl x /K? so = q zz against T in Fig. 3b. 
We find no appreciable change in q zz for 2-300K. This lack of 
change signifies that the localised character of the a\ g orbital 
persists to the Fermi liquid state across T*. Although q zz of 
the /-electron system has not been reported, the nonlinear re- 
lationship in the K~x plot may be a manifestation of Kondo 
coupling in the Ce and U-based compounds^"—. 

q zz reflects the 3d orbital polarisation when x a and x e sca l e 
to electron occupations. From q zz = — 0.7 we can evaluate the 
mixing ratio of the a\ g and e' g orbitals in L1V2O4. The singly 
occupied a\ g orbital has q zz = — 1.2, as observed in the insu- 
lating pyrochlore material LU2V2O7 (Ref. LL8|). In contrast, 
q zz vanishes for equivalent mixing of a\ g and e' g , as observed 
in a less correlated metal V2O3 (Fig. 3b). The observed in- 
termediate q zz in LiV2C>4 manifests a significant e' g contribu- 
tion to the spin susceptibility in throughout the temperature 
range. Namely, the e' g spin must be polarised via Hund's rule 
coupling to the localised a\ g spins under the magnetic field, 
although the itinerancy of e' g is much better than that of a\ g . 
The occupation ratio is evaluated as a\ g : e' g = 4 : 1 (see Meth- 
ods), corresponding to the electron numbers of n ai ~ 1 and 
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FIG. 3: Temperature dependence of spin and orbital polarisa- 
tions in LiV2C>4. (a) Local spin susceptibilities probed by 51 V 
isotropic shift K? SQ (circles with error bars), axially anisotropic shift 
— K^ x (triangles with error bars), and 7 Li Knight shift K s (green 
dotted line), where error bars are defined by the standard deviation. 
The right-hand axis shows the bulk susceptibility % (solid curve). A 
constant Kq = 0.4% obtained from the K-% plot (Fig. SI) is sub- 
tracted from K. (b) The effective 3d orbital polarisation q zz (blue 
circle with error bars) in comparison with the fully a\ g polarised 
value q zz £ = — 1.2 in LU2V2O7 (black square) 18 , and the fully un- 
polarised value q zz p = in a weakly correlated metal V2O3 (black 
circle) (Y.S., M.I., & Y. Ueda, unpublished results). q zz is obtained 
from K^ x divided by K^ o to cancel the spin polarisation. 



n e i ~ 0.25 for 3d 15 . The half-filling a\ g occupation is dis- 
tinct from that expected in the tight-binding calculation with- 
out electron correlations, where n aig : = 1 : 4 (Ref. Il2l) . 
However, it is compatible with the strongly localised a\ g pic- 
ture due to the strong renormalization into the Mott insulat- 
ing state 5,9,10 and provides microscopic evidence for orbital- 
dependent localisation, which is robust across T*. 

Dynamical spin susceptibility of the orbital-selective 
spin liquid Another interesting issue is the dynami- 
cal part probed by the nuclear spin-lattice relaxation 



rate T, 



(T\T) 1 is generally given by (T\T) 



£^|A ± (q)| 2 Im^ ± (q,co)/(0 (Ref. |19|), where Aj_(q) is the 
wave vector q component of the hyperfine coupling constant 
normal to the quantization axis, and ^(q, (o) is the transverse 
dynamical spin susceptibility at the NMR frequency (0. In a 
cubic lattice, A± (q) and x± (q? G>) are isotropic for the 51 V and 
7 Li sites. (T\T)~ l measured for 51 V and 7 Li (Fig. 4a) follows 

the linear relationship ( 51 7i!T) 1 = C( 7 T\T)~ l + Co, where 
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FIG. 4: The dynamical spin susceptibility compared to a static 
one in LiV204. a, The temperature dependence of the nuclear spin- 
lattice relaxation rate divided by temperature T, (T\T)~ l , obtained 
for 7 Li (solid circle with error bars) and 51 V (open square with er- 
ror bars) NMR (left axes), and the isotropic 51 V Knight shift 
(solid line, right axis) in L1V2O4, where error bars are defined by 
the standard deviation, b, C( 7 7ir) _1 (solid square with error bars) 
divided by K* SQ , where C( 7 7\r) _1 is equivalent to the spin part of 

( 51 t x t)- x . 



the linear coefficient C = 1.0 x 10 3 is close to the square ra- 

l.lxl0 3 andC = 



/51^1 A s \2 

tio of the hyperfine coupling, (1 7 , S 1S( ^ 



64 s _1 K arises from the ^-invariant orbital component. The 
scaling relation indicates that the Li sites probe spin fluctu- 
ations via the transferred hyperfine interaction and allows us 
to evaluate unobservable ( 51 7i!T) -1 at low temperatures from 

In the present case, (T\T)~ l is governed by paramagnetic 
fluctuations of local moments at high temperatures. Above 
150 K the scaling behaviour between ( 7 7ir) _1 and K? so is 
indeed observed (C( 7 TiTK? so )~ l = constant in Fig. 4b). Be- 
low 150 K a progressive ( 7 7i!T) -1 increase indicates antifer- 
romagnetic correlation, consistent with the growth of #(q, (0) 
at a finite q (= 0.64 -1 ) in the inelastic neutron scattering mea- 
surements below 80 K (Ref. [16|). Therefore the suppression of 
X s at low temperatures likely comes from the short-range anti- 
ferromagnetic correlation with the exchange interaction Jaf ~ 
150 K. Nevertheless, no long-range magnetic ordering occurs 
down to 1.5 K, the energy scale of ~ /af/100- It suggests 
that the frustrated a\ g spins form in a quantum liquid at low 
temperatures with low-lying excitations. 
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Discussion 

Our results provide significant insights into the formation 
of 3d HQ. As observed in the ^-independent q zz , the a\ g spins 
likely remain incoherent, even entering into a coherent 'Fermi 
liquid' state, and couple ferromagnetically to e' g spins. No 
indication of Kondo coupling was observed down to low tem- 
peratures despite the large antiferromagnetic fluctuations. The 
remaining local moments can be highly frustrated and carry 
large residual entropy 6-8 . The itinerant e' g electrons interact 
with the underlying spin liquid via the Hund's rule coupling. 
Thus the 3d HQ behaviour in L1V2O4 could be mapped on the 
frustrated ferromagnetic Kondo lattice. 

In the absence of antiferromagnetic Kondo coupling, the 
HQ formation has not established theoretically. In a Hubbard 
model calculation, the Kondo-like coherence peak appears on 
the boundary of the orbital- selective Mott transition for the a\ g 
part 10 . Even in such a case, % a may vary across T*, while % e 
be invariant. Our results suggest that, if the Kondo-like peak 
appears, a large fraction of the incoherent spins still remains 
and carries entropy. Such fractionalisation of the nearly lo- 
calised electron might be common to strongly correlated elec- 
tron systems 23 where localised and itinerant characters coex- 
ist. Furthermore, the chirality degrees of freedom might pro- 
vide appreciable entropy in the pyrochlore lattice 24 and con- 
tribute to the anomalous Hall effect 3 . 

The ORNMR technique offers new holographic experi- 
ments that could give microscopic insights into strongly cor- 
related electrons. Various orbital-resolved tools, such as X- 
ray absorption and photoemission spectroscopy, have been re- 
cently developed. The Knight shift measurement has a unique 
advantage in detecting the orbital-dependent local spin sus- 
ceptibility via the magnetic hyperfine interactions between d 
spins and on-site nuclear spins. The method has not been 
achieved in rare-earth heavy fermion compounds because of 
the difficulty in detecting NMR signals for the on-site nuclear 
spins 20 . Additional technical improvements in the NMR mea- 



surements may reveal the hidden orbital- selective Mott tran- 
sition in transition metal oxides, such as ruthenates, pnictides, 
and manganites. 



Methods 

NMR measurements. The ORNMR experiments were 
performed on a single crystal of L1V2O4 synthesised by the 
self-flux method 25 . The crystal with the octahedral shape was 
placed on a two-axis goniometer and rotated under a fixed 
magnetic field H = 9.402 and 8.490 T. The NMR spectra 
were obtained from spin-echo signals after two n/2 pulses 
separated by a time T. The 51 V NMR measurements were 
made only for powder samples above 50 K (Ref. [Tsb likely 
due to the fast spin-echo decay time T2 at low temperatures. 
To overcome this problem, we used a short T = 3 — 10/ls and 
a magnetic field precisely (<0.1°) aligned to the crystal axis 
equivalent to the magic angle of the nuclear quadrupole inter- 
action. Otherwise, the NMR signals were depressed owing to 
the fast 72. 

The angular dependence of the 51 V Knight shift K(6) with 
the local trigonal symmetry is fitted into the general formula 26 

K = K iso -&^sm26 (2) 
for the VI site and 

K = K iso ± fc^cos(20 ± 54.7°) (3) 
for the V2 or V3 site related by a mirror symmetry, where 

v K Z +2K X 
*nso — 3 

Magnetic hyperfine interactions. Magnetic hyperfine in- 
teractions between 3d electron spins and nuclear spins, and 
the analysis of the 51 V Knight shift in L1V2O4. Magnetic hy- 
perfine interactions in 3d systems are generally given 



1/ • I — KSi I 



(2/,-l)(2/j+3) 



{ \ (h • St) (h • I) + § (h • I) ft -Si)- l a (h + 1) (I • Si) } 



(4) 



where 1 ; , S/, and I denote operators of orbital and spin of the 
z-th electron, and the nuclear spin, respectively, the coefficient 
= 2^y n h(r~ 3 ) using the Bohr magneton /Ib, the nuclear 
gyromagnetic ratio y n , the Plank's constant h, and a radial ex- 
pectation value (r -3 ). The first term represents the orbital 
contribution that quenches in crystals but partly revives via 
the Van-Vleck process under the magnetic field. The second 
term arises from a Fermi contact interaction due to the core 
polarisation of inner s spins, giving the isotropic hyperfine 



coupling constant A? 



-&K8ij, where k is a dimension- 



less parameter (k ~ 0.5 for vanadates 18 ). The third term de- 
notes anisotropic dipole interactions determined by 3d orbital 



occupations, where the principal components are expressed as 
= -^^(—(lap +hfs! a p)/l by using the equivalent opera- 
tor of 3d angular momentum, 



<7a/3 = ^(W)3+^^a)-T/(/ + l)5 ai 3, 



(5) 



with the spin-orbit coupling parameter A, and the second- 
order matrix elements between the ground and excited states, 
. In the LS-coupling the sum of the terms for several elec- 
trons can be replaced by 



4 = ^LJ-^S-I-^ 



-§{(L-S)(L-I) + (L-I)(L-S)}-§L(L+1)(I-S) 



(6) 



where L and S are total orbital and spin, respectively, and £ = ^ 2/ _ p/+3) S (2L- 1 ) = IX ^ or 
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The experimental observable is the Knight shift tensor K = 
(K x ,K y ,K z ) defined as the resonance frequency shift due to 
the hyperfine interaction of Eq. (4). The spin component ^ s 
is obtained by subtracting the small orbital component Kq in- 
cluding the chemical shift and the Van-Vleck shift from the 
K-% plot in Fig. SI. For a paramagnetic system, S; in Eq. 
(4) is replaced by the effective electron spin polarization S; 
proportional to spin susceptibility X s - ln a multi-orbital sys- 
tem, X s is composed of orbital-dependent spin susceptibilities. 
Then K s is expressed by the sum of hyperfine fields from 3d 
spins. Whereas the isotropic part of K\ K* so (2K* + K*) /3, 
is given by 



^iso 



Nfl B 



the anisotropic part, K^ x = 2(K* — K*)/3, is expressed as the 
arithmetic average of the orbital-dependent spin susceptibil- 
ity weighted by the principal z component of the hyperfine 
coupling tensor, as shown in Eq. (1). In LiV 2 4 with a lo- 
cal trigonal distortion, 1.5 electrons are filled in two orbitals, 
a\ g and e f g , whose q a p are equivalent to those of d 3z i_ r i and 
d x 2_ y 2, respectively, taking the principal axes along the trigo- 
nal axis: q xx = q yy = 1, q zz = —2 for (^lg) 1 , while the values 
are numerically the same but reversed in sign for (e' g ) 2 . q aa 
vanishes when the two orbitals are equally occupied. Using 



the relation A a z 



, K* x is expressed as 



Kl 



(8) 



for the negligible AAj - as expected from the small Kq. To 
experimentally obtain the effective orbital polarization q zz , we 
take a ratio of Eqs. (7) and (9) and cancel out the numerical 
constants and X s - Namely, 



^^~(2/-l)^fe. 



(10) 



Here q zz = —8/7 for a fully a\ g polarized case (/ = 1), close 
to the experimentally obtained q zz ~ —1.2 in LUV2O7 (Ref. 
18). From the experimental result q zz ~ —0.7 in L1V2O4, we 
obtained / ~ 0.8, corresponding to the occupation ratio a\ g : 



C 7 ) e g = 4: 1 and hence n ~ 1 for a\ g and n ~ 0.25 for e' g . 



Electric hyperfine interactions The electrostatic hyperfine 
interaction can be a direct probe for 3 d orbital order. In the 
presence of the anisotropic electric field gradient around the 
nuclear spin, the 51 V NMR spectrum is split into seven lines 
for / = 7/2. The NMR spectrum becomes sharpest at [001], 
identifying the magic angle where the nuclear quadrupole 
splitting frequency <5v vanishes. Then <5v should have a max- 
imum at 60 = 54.7° satisfying <5v ~ (3cos 2 #o — 1) = 0, which 
exactly equivalent to the local trigonal symmetry. We ob- 
served a quadrupole splitting frequency Sv x = 90 kHz in the 
51 V NMR spectra at Ho || [110], by using a very short pulse 
interval time 1 = 3 /is. From the lattice symmetry, we can 
obtain Vq = 28 v x = 180 kHz. We confirmed that §v x was in- 
dependent of temperature down to 2 K and hence the orbital 
occupation was invariant across r*. 



The good linearity in K~x plots indicates X s = % a + % e ~ 
fX s + (1 ~f)X s > where / is the fraction of x a m X s - Tnen 
£L can be further reduced to 



^X s (2f-l). 



(9) 
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